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Integral form of the space-time-fractional Schrodinger equation for the scattering problem in the fractional 
quantum mechanics is studied in this paper. We define the fractional Green's function for the space-time 
fractional Schrodinger equation and express it in terms of Fox's H-function and in a computable series form. 
The asymptotic formula of the Green's function for large argument is also obtained, and applied to study the 
fractional quantum scattering problem. We get the approximate scattering wave function with correction of 
every order. 

1 INTRODUCTION 

The fractional calculus fl]|2], which is a generalization to the standard (integer) one, has been successfully applied in 
many science and engineering fields, such as anomalous transport, viscoelastic material, signal analysis and processing, 
control system |3}j6]. In recent years, the fractional calculus enters the world of quantum mechanics. In the fractional 
quantum mechanics |7, 8|, the master equation is the fractional Schrodinger equation [9 | instead of the standard one. The 
standard Schrodinger equation II 1 01 1 1 11 was reformulated by Feynman and Hibbs [12| using the path integral approach 
considering the Gaussian probability distribution. The Levy stochastic process is a natural generalization of the Gaussian 
process. The possibility of developing the path integral over the paths of the Levy motion was discussed by Kac lfl3ll . 
who pointed out that the Levy path integral generates the functional measure in the space of left (or right) continued 
functions having only discontinuities of the first kind. Then, Laskin Q generalized Feynman path integral to Levy one, and 
developed a space-fractional Schrodinger equation containing the Riesz fractional derivative ||2][9]. Then, he constructed 
the fractional quantum mechanics and showed some properties of the space fractional quantum system [8,9 14, 15 1. 

Afterwards, Naber lTT6l constructed a time-fractional Schrodinger equation by introducing the Caputo fractional deriva- 
tive 0] instead of the first order derivative over time to the standard Schrodinger equation to describe non-Markovian evo- 
lution in quantum physics. The Hamiltonian for the time fractional quantum system was found to be non-Hermitian and 
not local in time. Naber solved the time fractional Schrodinger equation for a free particle and for a potential well. Proba- 
bility and the resulting energy levels are found to increase over time to a limiting value depending on the order of the time 
derivative. More recently, from the standard Schrodinger equation, Wang and Xu fT7l established a Schrodinger equation 
with both space and time fractional derivatives, and solved the generalized Schrodinger equation for a free particle and for 
an infinite rectangular potential well. Then, a similar space-time-fractional Schrodinger equation is obtained by Dong and 
Xu iTTSTl from the space-fractional Schrodinger equation. They expressed this fractional equation in a more simple form, 
and studied the time evolution behaviors of the space-time-fractional quantum system in the time-independent potential 
fields. 

At present, the progresses for the space-time-fractional quantum system is fewer. Besides the results given in Refs. [ 17 , 
18] mentioned before, Jiang [20 1 developed a time-space fractional Schrodinger equation containing a nonlocal term, and 
obtained the time dependent solutions in terms of the H-function. All of these results are in the one-dimensional case. 
This paper focuses on the time-dependent scattering problem in the fractional quantum system described by the space- 
time-fractional Schrodinger equation, given by Dong and Xu, in the three-dimensional case. We will define the Green's 
function for the time-dependent 3D space-time-fractional Schrodinger equation for the scattering problem, and the Green's 
function will be calculated in terms of Fox's H-function and in a computable series form. The asymptotic property of the 
Green's function will also be discussed and applied to the fractional scattering problem. 
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2 GREEN'S FUNCTION OF THE FRACTIONAL SCHRODINGER EQUA- 
TION 

The space-time fractional Schrodinger equation [9 1 obtained by Laskin reads (in two dimensions) 

{iKftiP, tf/(r, t ) = M'cdix, t), ( 1 ) 

where i/f(r, f) is the time-dependent wave function , and 

= M [-D a (hVf + V(r, t)]. (2) 

Here, 

HP 

M - 



^p 1 p 



, D a with physical dimension [D a ] = [Energy] 1 a x [Length]" x [Time] " is dependent on a [we have D a = l/2m for 
a = 2, m denotes the mass of a particle] and (hV) a is the quantum Riesz fractional operator []2][7] defined by 

(W»(r,?) = " — t f d 3 pe*' r/ *|p| a f e- ,pr/ V(r, f)d 3 r. (3) 

(2tt#) 2 J J 

Note that by use of the method of dimensional analysis we have given a specific expression of D a in Ref. Ifl8ll as D a = 
c 2 ~ a /(am a ~ l ), where c denotes the characteristic velocity of the non-relativistic quantum system. 
Now, we define a Green's function of the FSE by 

[{ihflf, + £)(ftVf]G(r, f; r', t') = 6(r - r')6(t - t'), (4) 

with the causality condition 

G(r,f;r',f') = 0, when r<r'. (5) 
Here, T) — M ■ D a . Then, the space-FSE (Q} becomes an integral equation, 



i//(r,t) =ifr (r,t) + D jj G(r,t;r' ,t')V(r' ,t')i//(r' ,t')d 3 r'dt' , (6) 

in which 0"o(r, f) satisfies the free-particle Schrodinger equation, 

d 

[ffi-+B(fiVf]«r,^0. (7) 
at 

By use of the method of separation of variables, the basic solution to Eq. (0 can be easily obtained. 

^o(r, = e' (k ' r ~ a)//i , (a constant product factor is omitted ) (8) 
where E denotes the energy of the free particle, and k = (k x , k y ), in which k x , k y are arbitrary constants but satisfying 
|k| = Ah\ + = (E/D) l ^ a . Replacing k by momentum p, and E by D\p\ a respectively, the free particle solution (Ao(r, t) 
can be changed to the fractional plane wave solution [ 14 1, 

(Ao(r,f) = e' (pr - 2)|pl "' )/s . (9) 

Now we turn back to solve Eq. ©. Defining the Fourier-Laplace transform pair, with respect to r and f, of the Green's 
function G(r, f; r', t') as 

d 3 r J dfe" ,pr/ ^'G(r,f;r',f'), (10) 

^le'P^ + "G(p,.;r',f'). (11) 
(iTihy J c _ ioo 2m 
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After taking the above Fourier-Laplace transform with respect to r and f, Eq. (|4]i can be changed to 

(ihf^Gip, s; r', ?') + / _1 G(p, 0; r', f')] - D\p\ a G(p, co; r', ?') = e _fp * ' lh ~ s< , (12) 

in which G(p, 0;r',f') is the Fourier transform of G(r, 0;r',f') with respect to r. Then, taking account of the causality 
condition given by (0, we have G(r, 0; r', f') = 0, so that G(p, 0; r', f') = 0. Now, we can obtain 

g-ip-r'/fig-sr' 

G(p,w;r',f') = -——J — — — . (13) 
F {iKfsP - £>|p| a 

Inverting the Fourier-Laplace transform gives 

J_P efp-fr-n/n — — 1 . (14) 

To calculate the integrals in the above formula, more work is needed. Firstly, the s integration can be calculated by use of 
the residue theory, we have 

— I - i s = Res I - s sA= — = - M-otf-V/fi (15) 

2m (ihfsP-DW \(ihfsP-!DW ' °J pahf^ 1 iti/m^W ' K > 

Here, Res {*, s, so) denotes the residue [22] of *, with respect to s, at the unique pole so, and so — (flpD'^/O^)- Now 
Eq. (TBI) becomes 

G(r,?;r',?') = fe^^V^Wlpl^d'p, t>t', (16) 

(27rft) J J 

in which 

Nl = WH^-W ? = O 1 ^ - 0/«. v = «/A 7 = 0.08- 1)//J. 

To execute the above integration, we choose the spherical coordinates (p, 6, (p), with the positive direction of the p-axis 
along r - r'. Then, p ■ (r - r') = p\r - r'| cos 6, in which p and |r - r'| denote the magnitudes of the vectors p and r - r', 
respectively. Thus, Eq. ( [ToT l is converted into 



G(r,f;r',f') = l — \ d<9 dip \ e''' j|r - r ' |cosfl/ V'' flpl> |p| 2_y sinfld/? 

{2nfry Jo Jo Jo 

N 

2n 2 h 2 \r-r'\ J 



(17) 



In the last integral, substituting p for ^' v p gives 

G(r,t;r',t')=N--^-, (18) 
|r-r'| 

in which 

p l ~v sm{px)er if dp, x>0, (19) 

o 

|r - r'l Ni^ y - 2)/v 

Using the Mellin transform 11231 and its inverse transform, I(x) can be expressed in terms of Fox's H-function. Taking 
Mellin transform to I(x) with respect to x yields 



X+oo / j~+oo \ 

p x - y e- ip¥ I J %m{px)x s - 1 dx\ dp 



(21) 



= F(s)sin(y) | p x - y - s Q- ipV dp, (Res<l). 
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Then using the contour in Fig. Q] the above integral can be calculated, 



Thus, 



where 



Note that the formulas J23) 



f 

Jo 



^-r-^-'P'dp = -r(— 7 3 



-)e~^> 



I(s) = -r(sW( 2 7 V nffie- 2 ^" = -Ui(s)-il 2 (s)l 
v v 2 v 



2 — y — s ns 2 — y — s 

h(s) = r(s)T( r - )sin(-)cos(— f n) 

v 2 2v 



*- 2 ro)r( 



2—y—s 



r(f)r(i-f)r(i + ^)r(i-^£)' 

2 — y—s ns 2 — y—s 
his) = T(s)T{ y - )sin(-)sin(— f n) 

v 2 2v 



;r 2 r(s)r(^) 



r(f )m - |)r(^)r(i 



2-y-s \ ' 
2v > 



(22) 



(23) 



(24) 



(25) 



sin(7rz) = 



nz)T(i- Z ) 



, and cos(7rz) = 



r(i/2 + z)r(i/2-z) 



have been used here. Inverting the Mellin transform and comparing the expression with the definition of Fox's //-function 
(see Eqs. (lAlb and (IA2I) in the Appendix of this paper or Refs. [ 24-26 1), we obtain 



i r +; 

2m Jc-ic. 



v 



H 



3.3 



iH\j 



I(s) X - s ds 

(1 - (2 - r )/v, l/v), (0, 1/2), (1/2 - (2 - y)/(2v), l/(2v)) 
(0, 1), (0, 1 12), (1/2 - (2 - r )/(2v), 1 /(2v)) 
(1 - (2 - y)/v, l/v), (0, 1/2), (1 - (2 - y)/(2v), l/(2v)) 
(0,l),(0,l/2),(l-(2-y)/(2v),l/(2v)) 



(26) 



=/r 2 [Hx(x v )-iH 2 (x v )\, 
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where 



H 2 (x) = H 



(1 - (2 - y)/v, 1), (0, v/2), (1/2 - (2 - y)/(2v), 1/2) 
(0,v),(0,v/2),(l/2-(2-y)/(2v),l/2) 
(1 - (2 - y)/v, 1), (0, v/2), (1 - (2 - y)/(2v), 1/2) 
(0,v),(0,v/2),(l-(2- r )/(2v),l/2) 



(27) 
(28) 



Note that the first property of the //-function in the Appendix has been used (or see Property 1 .4 on page 12 of Ref. [24] ). 
Finally, we get 



G(r , f;rV) = ^l. 

|r-r'| 

£)(y+v-2)/(jSv)- 1 £(2-y)/v-3 

~2/?(f-?') (2 ~ y)/v |r-r'|/ 



Hi 



:r-r-| 



iH 2 



&iHt - t')h v 



ir-r'l 
- iH 2 



(29) 



t > f . 



3 COMPUTABLE SERIES FORM OF THE GREEN'S FUNCTION 



Using the series form of the Fox's //-function (see Property 3 in the Appendix of this paper or §3.7 in Ref. 12511 ). a 
computable form of the Green's Function can be obtained. For the //-function H\(x) given by Eq. (1271 1. we have 



r((2-y)/v + fc/v) 



(-1)* x k > v 



r(l + fc/2)r(-£/2)r(l/2 + (2 - y)/(2v) + fc/(2v))r(l/2 - (2 - y)/(2v) - k/(2v)) k\ v 



1 CO 

Ayr 



2 - y + k\ sin(-kn/2) cos[(l/v+ (k - y)/(2v))?r] (-1)* x 1 



k ^kjv 



k\ V 



2n + 3 - y 



cos 



2n + 3 - y 
7 

2v 



(-D" 
(2«+ 1)!" 



r (2«+l)/v 



In a similar way, we obtain 



1 ^ /2n + 3-y\ 

n=0 N ' 



sin 



2« + 3 - y 



(-D" 
(2« + 1)! 



(2n+\)lv 



Then, the Green's Function given by Eq. (1291 can be expressed in a series form as 



. . JV ^ /2n + 3-y\ / 2n + 3 - y \ (-1)" / |r - r'|\ 2 " 



(30) 



(31) 



(32) 



When j6 = 1, this formula reduces to the space-fractional case, see Eq.() in Ref. . 



4 ASYMPTOTIC PROPERTY OF THE GREEN'S FUNCTION AND ITS AP- 
PLICATION TO THE SCATTERING PROBLEM 

In quantum mechanics, the scattering theory H10II111 is often used to study the inner structure of a matter. The frac- 
tional quantum mechanics is a natural generalization to the standard quantum mechanics, so the research on the generalized 
quantum scattering problem under the framework of fractional quantum mechanics is meaningful. In the scattering prob- 
lems, we usually consider the behavior of the particles far away from the scattering center, and assume the potential V(r) 
is non-zero only in a small domain. Therefore, for r and r' in Eq. (0, we have [r| >> |r'|, or |r - r'| — > oo. In this section, 
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we study the asymptotic properties of the Green's function when |r - r'| — > oo. In this case, x = |r - r'\(^fi) l , defined by 
Eq. ( 1201 ). also approaches infinity. Recalling Eq. 129) . G(r, t; r', f ) can be rewritten in terms of x as 



Ntt 



G(r, t; r', f) = - ■ [77, (x v ) - iH 2 (x v )] , 



(33) 



in which the two H-functions H\(x) and H\(x) have been defined by Eqs. d27b and d28T i respectively. It is easy to verify 
that the two H-functions H\{x) and H\{x) satisfy the conditions needed by Property 4 in the Appendix. Thus, for large x, 
using the asymptotic formula dAlOt , we obtain 



M 



-r(2 - y) sin(y7r/2)x^~ + — cos[- +(v- 1)(— ) 1/(v_1) ]x^ +o(x^), 

7T Z HI 4 V v \ / 

1 rz? M 
#2(*) = - 7T(2 - y) sin(y7r/2)x » : sin| 



7T- 



(34) 
(35) 



in which 



Then, from Eq. 



„(2y-3)/[2(v-l)] 



M : 



V2tt(v- 1) 

we get the asymptotic formula of G(r, ?; r', f') for large |r - r'|. 



G(r,f; r',f') = 



2f +v h 3+v i 



exp 



/(a - 1) 



a/(a-l) 



+ O 



(36) 



When a — 2, j3 — 1, after some calculations, Eq. 



G(r,f;r',f') = - 



reduces to 

3/2 



vh \liMt- t')i 



im\r - r 
exp | -7-7- — I + o 



2h(t - ?') 



' V2m|r - r'| 



V(f - t')h ) 



(37) 



which accords with the exact result in the standard quantum mechanics 0271 . 

In the scattering problem, we can use merely the first term of the asymptotic formula d36*i > for G(r, f; r', f'), then an 
approximate wave function for the scattering problem can be obtained. Let's invoke the Born approximation lfTT1l28l : 
Suppose the incoming plane wave is not substantially altered by the potential. Then, in Eq. (O, it makes sense to use 



iA(r',f')~iAo(r',f') = e'' (kr '- £ '' )/fi . 



(38) 



Then Eq. (O becomes 



\fj{r, t) « i// (r, f) + 



2ft 3 

i(k-r' - Et')/h 



' 4 ff \r-r'r l 
-iJJ ~l^ eXP 

V(r',f')dVdf'. 



j(a-l) 



|r - r' 



(39) 



Furthermore, considering |r - r'| « |r| = r, we can get 

Ae m/4 r v-l 



t//(r, t) ^ i// (r, f) + 



2ft 3+ 



II e +v 



exp 



a/(a-l) 



Ka- Dl^l + '(k ■ r' - 



V(r',f')dVdf'. 



(40) 



The second term of the right side of the above formula gives the approximate scattering wave function. 

We can also generate a series of higher-order corrections to the approximate wave function. From Eq. (0, we can build 
an iteration scheme for the wave function as 



^ w (r, t) = (/r (r, f) + jj G(r, t; r', t')V(r', r')^ (n_1) (r', f')d 3 r'df'. 



(41) 
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Note that 4> (0 \r, t) = if/o(r, t) = e' (kr E< " , and </> (n) is the nth-order corrections to the wave function. For a given potential 
function, using Eq. ( |4TT >. the analytical approximate solutions of every order can be obtained. In a series form, we have 



= 00 + 



j GV<p a + jj GVGVfo + JJJ GVGVGV<p a + ■■■ . (42) 
In each integrand only the incident wave function (</>o) appears, together with more and more powers of GV. 



5 CONCLUSIONS 

In this paper, the three-dimensional space-time-fractional Schrodinger equation with time-dependent potential is stud- 
ied. We define the Green's function of the STFSE for the time-dependent scattering problem in the fractional quantum 
mechanics, and the STFSE is converted into an integral form. We give the mathematical expression of the Green's func- 
tion in terms of Fox's H-function and in a computable series form. The asymptotic property of the Green's function for 
|r - r'| — > oo (or |r| >> |r'|) was also given. In this case, the Green's function acts like an exponential function (see 
Eq. (136*1)). Using this result, We obtained the approximate scattering wave function for the time-dependent fractional quan- 
tum scattering problem (see Eq. (l40l). A series of higher-order corrections to the approximate wave function were also 
given in Eq. (1421 . These results are useful for the time-dependent scattering problem in the fractional quantum mechanics. 
All of these results contain those in the standard quantum mechanics as special cases. 
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APPENDIX: FOX'S //-FUNCTION AND SOME PROPERTIES 

The Fox's //-function [24 25 ] is defined by an integral of Mellin-Barnes type [29] as 



(auA l ),(a 2 ,A 2 ),...,(a p ,A p ) 
(b 1 ,B 1 ),(b 2 ,B 2 ),...,(b p ,B p ) 



where 



U%i + Bjs) UU rq - at - Ajs) 

nf =n+1 not + A tS ) u% m+1 m-bj- b jS ) 



(Al) 



(A2) 



The contour L runs from c - ioa to c + ioo separating the poles of F( 1 - a,- - A/s), (i — 1 , • ■ • , n) from those of T(bj + Bjs), 
(j = 1, • • • ,m). Here we present some properties of the //-function used in our paper. In order to give the results, the 
following definitions will be used, 



7=1 
P 



i= 1 i=«+ 1 

9 P 



./=! 



j=m+l 



(A3) 



p-q 



7=1 



7=1 



7=1 



/=I 



The following properties of the //-function can be found in Refs. 
Property 1: 
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1 Tjm.n 



!>■<■! 



(a p ,A p ) 
(b p , B p )) 



jjm.n 
H P,q 



(flp, kAp) 
(b P ,kB P ) 



. for k > 



Property 2: 



(j TTin,n 

Z P,1 



z 


(a p ,A p ) 


Tjin.n 
~ H p,q 


Z 


(a p +o-Ap,A p ) 


. for <x e 




(bp, Bp)) 






(bp + crB p ,B p ) 





Property 3: Explicit Power Series Expansion 

For A > , z + or A = 0, \z\ > 5, there holds the following expansion for the //-function [25], 



U]=ij*h T (bj - B jShk ) nil rq - q/ + AjS Hk ) (-i)* ^ 

where = (£>/, + k)/Bi u if the following conditions are satisfied: 

1. The poles of the gamma functions F(l - a, - A,s), (z - 1, ■ • • ,n) and those of Y{b j + Bjs), (j 
coincide: 

Ai(bj + l)± Bjiflj - k - 1), (i = 1, • • ■ , n; j = 1, ■ ■ ■ , m; k, I = 0, 1, 2, • • ■ ). 

2. The poles of the gamma functions T(bj + Bjs), (j - 1, ■ ■ • , m) are simple: 

+ * Bj(bi + k), (i + j; i, j = 1, • ■ ■ ,m;k, I = 0, 1, 2, • ■ • ). 

Property 4: Asymptotic Expansions at Infinity in the Case A > 0, A* = 

When the poles of the gamma functions T(l - a,- - A,s), (z = 1, • ■ ■ , n) are simple: 

Ai(\ - Aj + 1)* A/1 - fl; + k), (i + j; i, j = 1, ■ ■ ■ , n; k, I = 0, 1, 2, ■ ■ ■ ), 
and the condition (IA71 > is satisfied, the //-function has the following asymptotic expansion 



H ;*(z) = £ [h iZ «"- 1)/A > + o (z (a <- 1)/A >)] + A^ +I ^ A (c exp [(B + Cz i/A )i] ■ 

!=1 

db«ap[-(B + C z 1/A )i]) + o(z^ 1/2)/A ). 



where 



1 fl" i nbj - Bjjai - 1)1 Ad rq - fly + A/a,- - 1)1 A t ) 

A i U% n+l naj - Ajim - 1)/A ; ) U% m+l HI - ft; + Bjiat - 1)/A ; )' 



Ap /Ay 

2mA I <J 



B = 



(2// + 1)7T 
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A\l/A 



A = (27r)C-* +1) / 2 A-" f] A7 fl ' +I/2 f] 2# 



-1/2 



c = (2m)' 



m+n—p 



exp 



!=1 

/' 



7=1 



and t/o = ( — 2tt/)' 



m+n—p 



exp 



Z fl '-Z^ 

7=1 

m 

Z a '-Z^ 



\/=n+l 



7=1 / 
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